function [co] = transform_standard_rhombohedral(co)
%TRANSFORM_STANDARD_RHOMBOHEDRAL Summary of this function goes here
%   Detailed explanation goes here

a = co.LatticeVectors(1, :);
b = co.LatticeVectors(2, :);
c = co.LatticeVectors(3, :);

aH = a - b;
%bH = b - c;
cH = a + b + c;

aHLength = norm(aH);
cHLength = norm(cH);

aH = aHLength*[1 -sqrt(3) 0]/2;
bH = aHLength*[1 sqrt(3) 0]/2;
cH = cHLength*[0 0 1];

aR = 2*aH/3+bH/3+cH/3;
bR = -aH/3+bH/3+cH/3;
cR = -aH/3-2*bH/3+cH/3;


co.LatticeVectors = [aR; bR; cR];



end

